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SUMMARY * 

With the use of Kerman' s integrated momentum equation for the 
boundary layer and data on the wall-shearing stress and heat transfer 
in forced- convection flow, a calculation was carried out for the flow 
and heat transfer in the turbulent free-convection boundary layer on a 
vertical flat plate. The calculation is for a fluicl with a Prandtl 
number that is close to unity. 

< 

A formula was derived for the heat-transfer .coefficient that was 
in good agreement with experimental data in the range of Grashof numbers 
from 10^® to 10-^. Because of the good agreement between the theoreti- 
cal formula and the experimental data the formula may be used to obtain 
data for higher Grashof numbers. The calculation also yielded formulas 
for the maximum velocity in the boundary layer and for the boundary- 
layer thickness. 


INTRODUCTION 

Recent developments in the field of gas turbines have revealed the 
need for data on heat transfer in turbulent free-convection flow at 
very high Grashof numbers (10^ to 10^ ) . For instance, in using the 
method of free-convection cooling for turbine blades, the centrifugal 
forces generate a free-convection flow of the cooling fluid in the 
blade passages ’that is within the preceding range of Grashof numbers. 
Free-convection flow is also superimposed on the flow of the cooling 
air in the hollow blades of air-cooled gas turbines and may consider- 
ably influence the heat transfer under certain conditions. The radial 
flow present in the cooled boundary layers on the outside of cooled 
turbine blades is also of a type similar to the free-convection flow 
in the blade coolant passages. 

In order to understand and evaluate such cooling processes, 
information on the heat transfer, on the boundary- layer thickness, and 
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on the velocities connected with free-convection flow is necessary. The 
knowledge of turbulent free-convection flow, however, is limited. 
Experiments on plane vertical surfaces give heat- transfer data up to 
Grashof times Prandtl numbers of 10-^ (summarized in references 1 and 
2). Griffiths and Davis (reference 3) determined in addition tempera- 
ture and velocity profiles. Watzinger and Johnson (reference 4) meas- 
ured heat transfer and temperature and velocity profiles in a super- 
imposed forced- and free-convection flow in a vertical tube. Brown and 
Marko (reference 5) show by theoretical considerations that a general 
relation exists between the Grashof number that characterizes free- 
convection flow and the Reynolds number that characterizes forced flow. 
Colburn and Hougen (reference 6) derive a formula for the heat transfer 
in turbulent .free-convection flow on a vertical plate. However, only 
the laminar sublayer in the immediate neighborhood of the wall was 
investigated in reference 6 and the thickness of this sublayer, made 
dimensionless by the shearing stress velocity and the kinematic vis- 
cosity, was assumed to be the same as in forced flow. This analysis, 
therefore, gives no information on the whole boundary- layer thickness 
and on the velocities. 

The problem of turbulent free-convection flow is approached herein 
using another method. Certain shapes are assumed for the temperature 
and velocity profiles in the free-convection boundary layer. In addi- 
tion, an empirical relation for the shearing stress on the wall and a 
heat-transfer coefficient derived from forced-convection flow are used 
to estimate the boundary- layer thickness, the maximum velocity within 
the boundary layer, and a heat-transfer coefficient for free-convection 
flow. No experimental free-convection measurements are used in the 
calculations. 


SHEARING STRESS AND HEAT FLOW IN F0RCED- 
. ' CONVECTION BOUNDARY LAYER 

The relations for forced-convection flow that are needed in the 
theoretical free-convection calculations are compiled in this section. 
It is known that for Reynolds numbers that are not too high, the 
velocity profile in a turbulent boundary layer on a flat plate can be 
represented by the equation 

1/7 
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(All symbols are defined in the appendix. ) The shearing stress on the 
plate surface in such a flow is given by the equation (reference 7), 


T 


W 


0.0225 



( 2 ) 


This equation is the relation for the shearing stress on the wall that 
is used in the derivation of the boundary- layer equations for the tur- 
bulent free- convection boundary layer. By introducing the expression 
for the velocity profile, the shearing stress can also be expressed by 
the velocity u-]_ at the outer edge of the boundary layer and by the 

boundary- layer thickness & 


0.0225 




(3) 


For fluids or gases with a Prandtl number equal to one, the tem- 
perature profile is similar in shape to the velocity profile. When the 
temperatures within the boundary layer are measured with the temperature 
outside the boundary layer as reference, the expression for the tem- 
perature profile is 

9 = *wt-(|) l/7 ] M 

Equation (4) is valid for values of y less than 6. For larger values 
of y, 0=0. 

Reynolds analogy between the turbulent exchange of momentum and 
heat gives the relation (reference 8) 


3 . = 

T 


dt 

gc p du 


(5) 


Because the temperature and velocity profiles are similar in shape for 
a Prandtl number equal to one the finite temperature and velocity dif- 
ferences at two arbitrary points in the flow can be introduced into the 
last equation instead of the differentials. Using the differences 
between the wall and the flow outside the boundary layer and specifying 
the heat flow and the shearing stress for the position at the wall 
give the following equation: 


w 
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6w 

gc P^ 


( 6 ) 
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The introduction of the law (equation (3)) for the shearing stress into 
this expression yields the following equation: 


q v — 0.0225 gpCpU^© w 



(7) 


Heat transfer is usually calculated with dimensionless moduli. Such a 
modulus is the Stanton number 


St 


% 

gpcpu 1 9 w 



Experimental investigations show that for fluids with Prandtl numbers 
from approximately 0.5 'to 50 (reference 1, p. 520) the . same relation 
holds when it is multiplied by the factor Pr - ^/^. 


St 


*w 


gP c p u l®w 



( 8 ) 


This relation for the heat flow q w is used in the derivation of the 
boundary- layer equations for the turbulent free- convection boundary 
layer . 


Replacing the boundary- layer thickness & by 



which is the boundary- layer thickness on a flat plate in turbulent 
forced-convection flow (reference 1, p. 481), transforms equation (8) 
into the widely accepted formula for heat transfer on a flat plate in 
the turbulent range 



O) 


DETERMINATION OF TURBULENT FREE- CONVECTION BOUNDARY LAYER 

Derivation of Boundary-Layer Equations 

If a stationary plane vertical wall is heated to a temperature 
higher than the surroundings, the layer of fluid adjacent to the wall 
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is heated by conduction from the wall. In this way, buoyancy forces 
are generated that cause this layer to flow in an upward direction. 

This layer of fluid adjacent to the wall to which the vertical motion 
is confined is called the free- convection boundary layer. The boundary 
layer begins with zero thickness at the lower end of the. vertical wall 
and increases in thickness in the upward direction. The boundary layer 
becomes turbulent, depending on the critical Grashof number, at a cer- • 
tain distance from the lower end of the wall. The distance measured 
vertically from the lower end of the wall is called x and the dis- 
tance normal to the wall y. In order to determine the boundary- 
layer thickness for steady state, a small stationary volume element 
in the turbulent region of the boundary layer is considered. Figure 1 
shows this volume element. The dimensions of the element are dx along 
the wall and 7 normal to the wall. The length l should be larger ' 
than the boundary- layer thickness 6. For two-dimensional flow, the 
dimension of the volume element normal to the plane of figure 1 may be' 
considered to be unity. The upward velocity of the fluid in plane 1-1 
at a distance y from the surface of the wall is u. Then the mass 
flow through a small area with a width dy is pu dy and the flow of . 
momentum in the x-direction is pu^dy. The momentum flow in the 
x-direction entering the volume element through plane 1-1 is 

J' 1 pu 2 dy 

In progressing to plane' 2-2 the momentum flow changes by 



The mass flow entering the volume element through plane 1-1 is 
generally different from the mass flow leaving the element through 
plane 2-2. Therefore fluid enters or leaves the volume element through 
the plane parallel to the wall at a distance l. Because it is assumed 
that the velocity in the x-direction outside the boundary layer is so 
small that it can be neglected, no momentum in the x-direction is 
carried through the plane. 

The rate of change of momentum must be in equilibrium with the 
forces acting in the x-direction on the fluid within or on the surface 
of the volume element considered. A shearing stress acts on the 

wall. The force connected with this stress is T w dx. No shearing 
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stress occurs on the surface of the volume element that is parallel to 
the wall and at a distance l from the wall because outside the bound- 
ary layer the velocity in the x-direction is zero. 

According to the boundary- layer theory, the pressure change can be 
neglected along any normal to the surface. A constant pressure differ- 
ence dp therefore- exists between planes 1-1 and 2-2. This pressure 
difference gives a force on the volume element of magnitude l dp. In 
addition to the forces on the surfaces there is a force due to the 
weight of the fluid within the volume element 



The process of summing up all the forces and equating them to the 
change in momentum flow gives the momentum equation - 



No flow exists outside the boundary layer and therefore the pres- 
sure difference between planes 1-1 and 2-2 is balanced by the. weight of 
the fluid layer between the planes 


dp = p 5 g dx 


When both sides of this equation are multiplied by l and the right- 
hand side of the equation is changed to. the integral form, the follow- 
ing equation is obtained 


l dp = 



P 6 g <*yJ 


dx 


Introducing the preceding equation into the momentum equation gives 


( 

p \ 


P 

1 

pu 2 dyj dx = g 


(p s -p) <Jy 

V 

0 / 


0 


dx - t dx 
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Introducing the expansion coefficient defined by the equation 




the first term on the right-hand side of the preceding equation can be 
transformed into 

» pz 

g dx Pp(t-tg) dy 

Jo 


Designating the difference between the temperature t at the distance y 
and the temperature tg outside the boundary layer by 6 , the following 
expression is obtained: 

Pi 


g dx 


Pp0 dy 


In the applications considered, the essential influence of the 
density changes on the flow is taken into account by the introduction 
of the expansion coefffcient (3. The other influences of variable density 
on the flow and the variation of the expansion coefficient ,6 with 
small temperature differences are negligible. Therefore, (3 and p 
can be assumed constant in the preceding expressions. The momentum 
equation for the free-convection boundary layer therefore 
becomes 


_d_ 

dx 



( 10 ) 


A similar equation is set up for the heat flow through the volume 
element in figure 1. The heat carried with the fluid through plane 1-1 
is 


gpCp I u0 dy 

| 

where the enthalpies Cp0 are measured from the temperature outside 
the boundary layer. The specific heat and density are considered 
constant. 
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The heat carried out of plane 2-2 differs from that carried into 
plane 1-1 by 

' "I 


gpc p as 


u0 dy / dx 


The difference in the heat flov through planes 1-1 and 2-2 must come 
from the surface and the heat flow leaving the plate per unit time and 
area is therefore, 



‘"PI 



(u) 


Equation (11) is the heat-flow equation for the free-convection boundary 
layer. 

Equations for the laminar free-convection boundary layer that are 
analogous to equations (10) and (ll) were derived and solved by Squire. 
(See reference 9.) 


Solution of Boundary-Layer Equations 

Equations (10) and (ll) are sufficient to calculate the houndary- 
layer thickness 5 and the velocity u when the shape of the velocity 
and temperature profiles within the boundary layer and the laws for ,the 
shearing stress and the heat flow on the wall are known. It was sug- 
gested by von Karman that an approximate solution may be obtained by 
introducing approximate shapes for both profiles (reference 7). The 
accuracy of the solution is better the more closely the assumed shapes 
correspond to the real profiles. The method, however, proved compar- 
atively insensitive to changes in the profile shape. 

As mentioned in the section " INTRODUCTION" some information on the 
shape of the velocity and temperature profiles in turbulent free- 
convection flow can be obtained from reference 3. In figure 2, the 
measured values of both profiles in the turbulent range are presented 
with the ratios 6/© w and u/u^ plotted against the distance from 

the wall in an arbitrary scale. The value u mny was determined by 

drawing curves through the measured velocity values given in table VI 
of reference 3. 
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The distance from the vail vhere the curves u/u mflT have the value 
0.5 was arbitrarily assigned the value 1. The lines in figure' 2 are the 
curves by which the temperature and velocity profiles are approximated 
herein. These profiles are presented in figure 2 in such a way that the 
velocity profile has the value 0.5 at a distance 1 and the temperature 
profile the value 0.2 at a distance 0.5 from the wall. An indication of 
the shape of turbulent free- convection flow profiles can also, be obtained 
from the measurements by Watzinger and Johnson (reference 4) on mixed 
forced- and free- convection flow in a vertical tube. Some of these 
measured profiles are presented in figure 3. 

Equation (4), which describes the forced- convection temperature 
profile, is used for- the temperature profile in turbulent free- convection 
flow. The profiles of figure 2 indicate that this assumption is rea- 
sonable. The velocity profile in free-convection flow, however, differs 
from the forced-flow profile by having a velocity of zero outside the 
boundary layer. Accordingly, equation (l) is multiplied by; a function 
of y/&, which brings the velocities back to the value zero at the 
outer edge of -the boundary layer where y = 6. Two simple expressions 
that fulfill this condition and that agree quite well with the meas- 
ured values are plotted in figure 2. The velocity profile that rep- 
resents the measured points somewhat better is used in the following 
calculations : 

» - (if 7 1 1 - 1)* <*) 

Equation (12) is valid for values of y less than 8. For larger 
values of y, u = 0. The maximum velocity u n of this profile can 
be found by differentiating equation (12), giving 


u max - 0.537 u-^ 


(13) 


By use of equations (4) and (12) and the fact that u = 0 and 0=0 
for distances greater than 5 from the wall, the integrals in the 
momentum and heat-flow equations (equations (10) and ( 11)) become 



dy = 0.0523 &u. 



dy = 0.125 80 w 
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r . ■ 

/ u0 dy = 0.0366 8^0^ 

JO 

The same value of 6 was used for the velocity and temperature 
profiles. Calculations on forced- convection heat transfer indicate 
that the same value of 6 for the velocity and temperature profiles 
can be used for fluids that have a Prandtl number close to one (ref- 
erence 8). For very large or small Prandtl numbers in forced flow, 
the thickness of the temperature boundary layer is considerably smaller 
or larger than the thickness of the velocity boundary layer. This con- 
dition is also probably true for free- convection flow. 

In order to solve boundary- layer equations (10) and (11) the laws 
for the shearing stress on the wall and the heat flov must be intro- 
duced. It is assumed that, in the layers very near the wall, the 
conditions are similar for free- and forced- convection flow and that 
the same laws for the shearing stress and the heat flow that are used 
in forced- convection flow can therefore be used in free-convection 
flow. A similar assumption for the shearing stress was used by von 
KArmAn to calculate the radial flow on a rotating disk, with satis- 
factory results (reference 7). 

As a consequence of the preceding assumption equations (2) and 
(8) will be used for free-convection flow. Very near the wall the 
second term in equation (12) for the free- convection- velocity profile 
is one and, therefore, the same transformation that led in forced flow 
from equation (2) to equation (3) can be made for free-convection flow. 
By introducing the evaluated integrals and equations (3) and (8) the 
momentum and heat-flow equations become 



Equations (14) and (15) are total differential equations from which the 
two unknown values 6 and U]_ may be determined as functions of x. 

Equations (14) and (15), can be solved by introducing 

U 1 = C u x 


(16) 
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6 = C 5 x n (17) 

i • 

The introduction of these values transforms equations (14) and (15) into 


0.0523(2nH-n) C u 2 C 8 x 2m+n-1 = O.125gp0 w C 6 x n - 


0.0225 C u 7/4 C 6 - 1 / 4 V 1 / 4 x (7m/4)-(n/4) 


0. 0366(m+n)C u CgX 


m+n-1 


= 0.0225 (Prf 2 / 3 S 1 / 4 C u 3/4 3m/4)-(n/4) 


Because equations (18) and (19) must he valid for any value of x, the 
exponents of x must he identical, that is, 


2m+n-l = n 


7m n 
4 ‘ 4 


, 3m n 

m.n-1 - T - j 

It can he seen that this condition is fulfilled if 

“ a »-n < 

With these values of m and n equation (19) can be solved for the 
constant c u 

C u =0.0689 VC 6 ' 5 (Pr)‘ 8 / 3 ( 

Introducing this constant into equation (18) gives 


C R iU = 0.00338 


[u0.494(Pr)2/ 3 ] (Pr)' 16 / 3 




Introducing the Grashof number Gr = gives for the velocity 

V 2 

u x = 1.185 2 (Gr) 1 / 2 [\+0. 494(Pr) 2 / 3 ] 
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and for the boundary- layer thickness 


6 = 0.565 x(Gr)" 1 / 10 (Pr)" 8 / 15 j[l+0.494(Pr) 2 / 3 ]] 


1/10 


(24) 


More significant than the value uq_ is the maximum velocity 
u max the boundary layer. (See equation (13).) A Reynolds 

number is obtained to make this velocity dimensionless 


Re max = "HT* '= °* 636 (G r ) 1//2 [l+0.494(Pr) 2 / 3 ] 


- 1/2 


(25) 


A more significant value characterizing the thickness of the bound- 
ary layer than the value S : used up to now is the displacement thick- 
ness 6*, which is defined for free -convection flow by the equation 

P>oo 


5* = 


d 


u 


0 %ax 


dy 


With the use of equations (12) and (13) it is found that 

6* = 0.272 6 

The displacement thickness made dimensionless by the distance x 
from the leading edge of the plate is therefore 


^ = 0.154 (Gr) _1 / 10 (Pr)' 8 / :iS [l+0.494(Pr) 2 / 3 ] 


l/lO 


The heat-transfer coefficient H = q w /© v at the point x on the 

wall is found by introducing the boundary- layer thickness S and the 
velocity u-^ into equation (8). Changing to the usually presented 
Nusselt number Nu = Hx/k gives 


Nu = 0.0295 (Gr) 2 ^ 5 (Pr) 7/ ^ 15 




1+0.494 (Pr) 


2/3J 


-2/5 


(26) 


In order to compare this equation with experimental results, it is nec- 
essary to change from the local heat- transfer coefficient to the average 
value along the plate. By introducing the expression for the Grashof 
number into equation (26), it can be seen that the local heat- transfer 
coefficient is proportional to the power 0.2 of the distance x 
(H = CgX^*^) . With the assumption that the boundary layer is turbulent 
from the leading edge, the average heat-transfer coefficient becomes 
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In reality the boundary layer is first laminar and only at a certain 
distance from the lower edge of the plate does it become turbulent. 

The preceding expression for the average heat- transfer coefficient can 
therefore be expected to represent the true value only at Grashof num- 
bers so high that the extent of the laminar boundary layer is small 
compared with the total length x. This limit for the Grashof number 
seems to be near 10^®. For Grashof numbers higher than this value the 
average Nusselt number can be calculated with equation (27); therefore. 


Nu 


= 0.0246 (Gr) 2 / 5 (Pr) 7 / 15 [l + 0.494 (Pr) 2 / 3 (28) 


av 


In order to determine to what extent the shape of the velocity 
profile influences the results, a calculation was made with the second 
velocity distribution shown in figure 2. 



(29) 


Only the numerical constants are influenced by the change in the profile 
shape- The constant preceding the value (Pr) 2 / 3 changes from 0.494 to 
0.342 in all equations. In addition, the value 0.636 changes to 0.487. 
in the equation for Be^x; in the equation for S*/x the value 0.154 
changes to 0.200 and in the equation for Nu av the constant 0.0246 
changes to 0.0198. 


COMPARISON WITH EXPERIMENTS 

In figure 4 the results of experiments carried out in different 
investigations (references 1 and 2) are plotted as the average Nusselt 
number Nu av against the product Gr Pr. For the lower values of 
Gr Pr the experimental results quite accurately fit the equation 

. Nu av = 0.555 (Gr Pr) 1 / 4 

which was, with a slight adjustment of the constant, theoretically 
derived for laminar free-convection flow. The experimental results in 
the turbulent range (Gr Pr = 10 10 to 10 12 ) may be represented by 
the equation 
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Wu aT = 0.0210 (Gr Pr) 2 / 5 (30) 

The exponent of the Grashof number in this equation is the same as that 
of equation (28) derived in the previous section. When equation (28) 
is transformed into the form of equation (30) in such a way that the 
values of both equations are the same for Pr = 0.72, the constant of 
equation (28) becomes 0.0210. The heat-transfer coefficient derived 
with equation (12) for the velocity profile is in perfect agreement 
with the experimental results. Such agreement is probably a coinci- 
dence. The profile given by equation (29), which does not fit the 
velocity distributions measured by Griffiths and Davis in reference 3 
as well as does the first profile, gives heat- transfer coefficients 
that are 17 -percent lower than the measured values. 

The values for the maximum velocity within the boundary layer and 
the boundary-layer-displacement thickness agree poorly with the 'values 
measured in reference 3 in the turbulent range. Heat- transfer coeffi- 
cients were not measured therein on the experimental apparatus on which 
velocity and temperature profiles were obtained. 

Whereas the velocity and temperature profile shapes as measured in 
reference 3 are typical of turbulent free-convection flow the order of 
magnitude of these profiles appears to be in error. It can be shown 
that there is disagreement within the measured values themselves. The 
heat given off by the plate to the air stream must be carried away 
Within the boundary layer. The measured temperature and velocity pro- 
files as well as the measured maximum velocity and the boundary- layer 
thickness vary little along the plate in the turbulent range , which 
means that only a small part of the heat given off by the wall is found 
in the boundary layer. The horizontal dimension of the plate may have 
been too small compared with the vertical length to make the flow in 
the center part two-dimensional and air may have flowed into the bound- 
ary layer from the sides. 


SUMMARY OF RESULTS 

With the use of Kerman 1 s approximate method, a calculation was 
carried out for the flow and heat transfer in the turbulent free- 
convection boundary layer on a vertical flat plate. The calculation 
used relations for the heat flow and shearing stress on the wall devel- 
oped for forced turbulent flow and velocity and temperature profiles 
that approximate well the shapes of profiles measured by Griffiths and 
Davis. 
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A formula was derived for the heat- transfer coefficient that was in 
good agreement with measured values in the range of Grashof numbers 
from 10^ to 10^ and that can be used to extrapolate the values into 
the range of higher Grashof numbers. The formula is valid for Prandtl 
numbers that are close to unity. 

The calculation also yielded formulas for the maximum velocity in 
the boundary layer and for the boundary- layer thickness. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for 'Aeronautics, 
Cleveland, Ohio, July 12, 1950. 
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APPENDIX - SYMBOLS 

The following symbols are used in this report: 
thermal diffusivity, sq ft/sec 

constant for variation of heat- transfer coefficient with x 0,2 

constant for variation of velocity u in boundary layer 
with x® 1 

\ 

constant for variation of boundary- layer thickness with x 11 
specific heat at constant pressure, Btu/(lb)(°F) 


Grashof number. 




acceleration due to gravity, ft/sec 2 
heat-transfer coefficient, Btu/(sq ft)(sec)(°F) 
heat conductivity, Btu/(ft) (sec) (°F) 
length (fig. 1), ft 
exponent 

Nusselt number, — 

’ k 


exponent 


gcnd 


Prandtl number, - = —E 
a k 

pressure, lb/(sq ft) 

specific heat flow, Btu/(sec)(sq ft) 


u x 
max 


^by Reynolds number based on maximum velocity ^ 

1/4 


St 


Stanton number. 


= 0.0225 (JL-) 

gpCpUi^w V u l 6 / 
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t temperature, °F 

u velocity component in x-direction, ft/sec 

velocity outside boundary layer of comparable forced- convection 
flow, ft/sec 

v specific volume, cu ft/lb 

x coordinate (distance along plate from starting point of bound- 

ary layer), ft 

y coordinate (distance from wall), ft 

3 expansion coefficient, l/°F 

8 boundary- layer thickness, ft 

8* displacement thickness of boundary layer, ft 

9 temperature difference, °F 

9 V temperature difference between wall and fluid outside of bound- 

ary layer, °F 

H absolute viscosity, pv> lb sec/sq ft 

V kinematic viscosity, sq ft/sec 

p mass density, (lb)(sec^)/ft’ 

T shearing stress , lb/sq ft 

Subscripts: 
av average value 

max maximum value 

w on wall 
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Figure 1. - Nomenclature in deriving equations for free -convection boundary layer. 
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Figure 3. - Temperature and velocity profiles in mixed forced- and free- convection flow of water in vertical 
tube of 50.3 millimeter (1.98 in.) diameter. Data from figure 10 of reference 4. 





Vertical cylinder l ( Ee ference 1, fig. 25-3) 

Vertical plane J 

9. 01- in. plate [ (Reference 2, fig. 129 
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Figure 4. - Average Nusselt number for free- convection flows on a vertical plate 
Experimental data from references 1 and 2. 



